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Abstract 



The perfect quenching of spin tunneling that has previously been discussed in 
terms of interfering instantons, and has recently been observed in the magnetic 
molecule Fes, is treated using a discrete phase integral (or Wentzel-Kramers- 
Brillouin) method. The simplest model Hamiltonian for the phenomenon 
leads to a Schrddinger equation that is a five-term recursion relation. This 
recursion relation is reflection-symmetric when the magnetic field applied to 
the molecule is along the hard magnetic axis. A completely general Herring 
formula for the tunnel splittings for all reflection-symmetric five-term recur- 
sion relations is obtained. Using connection formulas for a new type of turning 
point that may be described as lying "under the barrier" , and which underlies 
the oscillations in the splitting as a function of magnetic field, this Herring 
formula is transformed into two other formulas that express the splittings 
in terms of a small number of action and action-like integrals. These latter 
formulas appear to be generally valid, even for problems where the recursion 
contains more than five terms. The results for the model Hamiltonian are 
compared with experiment, numerics, previous instanton based approaches, 
and the limiting case of no magnetic field. 
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I. INTRODUCTION 



In several previous papers |], we have discussed the tunneling of a spin described by the 
model Hamiltonian 

U = -k 2 J 2 z + (h - k 2 )J 2 x - g/ 2 B J ■ H, (1.1) 

where J is dimensionless spin operator with components J x , J y , and J z , H is an external 
magnetic field, and k\ > k 2 > 0. Since most of the earlier work has been for the special 
case where H||x, let us first limit ourselves to that. Viewed as a function of a classical 
angular momentum vector J of fixed length J, this Hamiltonian has two degenerate minima 
for H x < H c = 2kiJ / g^B- On general grounds we expect quantum mechanical tunneling 
to lift the degeneracy, and split the energies. The surprise is that the ground state tunnel 
splitting, A , is an oscillatory function of H x , vanishing exactly whenever 
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where A = k 2 /ki, and n = 0, 1, . . . , 2J — 1. 

The quenching of A was first studied purely as a theoretical curiosity, and explained 
in terms of instantons |]l|(a). Since then, this effect has been observed in the magnetic 
molecule [(tacn) 6 Fe 8 2 (OH) 12 ] 8+ (or just Fe 8 for short), which is approximately described 
by the Hamiltonian flOP , with J = 10, and k x « 0.33 K, and k 2 » 0.22 K §-§]. Motivated 
by a desire to use only elementary methods of analysis, the problem was restudied 01(d) 
using the discrete phase integral (DPI) method (also known as the discrete Wentzel-Kramers- 
Brillouin method). An approximate version of this method was developed and applied to Fe§ 
independently by Villain and Fort || . The results of these later studies confirm Eq. ( |1.2j ) for 
the ground pair quenching points, and also find additional quenching points, as we discuss 
next. 

To help grasp the full richness of the spectrum of the Hamiltonian ( |1.1| ), we show in 
Fig. |l] the results of numerical calculation of the energies as a function of H x , for J = 3, 
for three different values of H z . In all three cases, H y = 0. In part (a), H z = 0, and we 
have the symmetric situation mentioned above. Note that (i) the lowest two energy level 
curves cross six times (including negative values of H x ), and (ii) the crossing points are 
perfectly periodically spaced, in complete accord with Eq. (|1.2|). Wernsdorfer and Sessoli 
have shown the existence of analogous crossings for Fe§. To quantitatively account for 
the observed period, one must include higher order anisotropics in the Hamiltonian. This 
does not change the basic physics. In addition, Fig. |l](a) also shows a number of crossings 
of higher energy levels, of which the analog in Fes has not been seen yet. 

In Fig. |l](b), H z has a specific non-zero value. The problem is no longer symmetric, and 
one of the classical minima is lower than the other. Correspondingly, we see that the lowest 
quantum mechanical state is always non degenerate. However, ignoring tunneling for the 
moment, the first excited state in the deeper well can have an energy equal to that of the 
lowest state in the shallower well if H z is correctly chosen. And indeed, we see from the 
figure that the second and third energy levels do cross a number of times. These crossings, 
when H has an easy component, were not anticipated in Ref. 0(a), and were discovered 
by Wernsdorfer and Sessoli experimentally. As seen in the experiments, the crossings in 
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Fig. |l](b) are shifted by half a period from those in Fig. |l|(a). Note that as in part (a), 
Fig. |l|(b) displays crossings between yet higher energy levels (the fourth and fifth, e.g.), 
which have also not been seen experimentally yet. 

This pattern continues as H z is increased still further [Fig. fj(c)]. Now the lowest two 
levels in the deeper well are nondegenerate, and the lowest crossings are between levels 3 
and 4. Comapred to Fig. 0(b), these crossings are shifted by yet another half-period, just 
as seen experimentally. Again, there are crossings between higher pairs of levels, and again 
only those between levels 3 and 4 have been seen in Fe§. 

It is clearly interesting to understand the structure in the energy spectrum analytically, 
and numerical diagonalization alone cannot provide this. When J is of order 10, as it is for 
Fes, a semiclassical analysis is natural, and it is profitable to think of the energy differences 
amongst low lying levels in terms of tunneling. Such an analysis was done in Ref. 0(a-d). In 
this paper, we shall elaborate on our earlier DPI analysis [0(d), and provide several results 
that are more generally applicable to Hamiltonians other than Eq. ( |1 . 1| ) . We shall limit 
ourselves, however, to problems which are analogous to symmetric double-well potentials in 



the continuum WKB case. In the context of Eq. ( |1 . 1| ) , this means that H||x. The cases where 
H y or H z are also nonzero correspond to an asymmetric potential, and will be considered in 
a second paper. 

Before describing the results of our analysis, however, let us digress to make two points. 
The first is the issue of degeneracy and its connection with symmetry in light of the von 
Neumann- Wigner theorem. When H||x, or H||z, Ti is invariant under a 180° rotation about 
x or z, so energies of levels that are odd and even under this operation can intersect. 
The quenchings for H||x [Fig. [l|(a)] can be understood as instances of this phenomenon 
[lj(b). When H has both x and z components, however, 7i has no symmetry, and the level 
crossings [Fig. [TJ(b), (c)] are nontrivial instances of conical intersections or diabolical points 
|7]||. Viewed in the larger H x -H z plane, or the full three-dimensional space of magnetic 
fields H, however, all points of degeneracy are diabolical. 

The second point is that there are several other special features in the spectrum, which 
are evident from numerical analysis for several different J, and can also be seen in Fig. |l|. 
First, the successive half-period shifts in the crossing fields as we go from (a) to (b) to (c) 
in Fig. m mean that the diabolical points form part of a centered rectangular lattice in the 
H x -H z plane. The length of the rectangular unit cell along H x can be read off Eq. ( |1.2|) , 
while that along H z is given by @,|] 

A 1 / 2 

AH Z = —H c , (1.3) 

where A = k%/k\. Second, at a diabolical point, we often find simultaneous degeneracy 
of more than one pair of levels to very high accuracy if not exactly. All these facts are 
captured by the DPI analysis. In fact, in the case H||x, all the available evidence to date - 
exact diagonalization for small J, perturbation theory in A = numerics — indicates 

that the simultaneity of the degeneracy of many pairs of levels, as well as the values of 
the degeneracy fields, are exactly given by the leading semiclassical analysis, i.e., Eq. ( |1.2|) 
0(e). These facts point to the existence of a higher dynamical symmetry, but that is not yet 
established. Further, when higher order anisotropy terms are included in the Hamiltonian 
to obtain quantitative agreement with experimentally observed period @, the numerical 
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evidence indicates that although the simultaneous degeneracy of several pairs of levels and 
the perfect lattice of diabolical points are no longer exact properties, they continue to hold 
to rather good approximation |K| . 

The plan of the paper is as follows. In Sec. [H] we outline the DPI approach. The 



basic idea is that in the J z basis, Schrodinger's equation for Eq. (|1 . 1|) has the form of a 



recursion relation or difference equation, as opposed to a differential equation for a massive 
particle in a one dimensional potential V(x). This difference equation can be solved in close 
analogy with the continuum WKB approximation. We will see that compared to previous 
DPI studies [TT|-|il| new types of turning points arise in the study of Eq. (|1.1|), because the 
recursion relation has five terms as opposed to three in the earlier studies. These turning 
points have no continuum analogue. Our present discussion will rely on physical arguments 



and correspondence with the continuum case. A more formal discussion is given in Ref. |15 



In Sec. [nj we develop an analogue of Herring's formula ||T6| , [T7| for problems leading to 
five term recursion relations. In the continuum case, for a symmetric double well potential 
[V(— x) = V(x)], this formula expresses the splitting for the nth pair of levels in terms 
of the x = values of the wavefunction ip n (x) and its derivative ip' n (x) for the nth state 
localized in one of the wells. In Sec. [IV] we will use the DPI method to find the analogous 
discrete wavefunction near the center of the potential, and use our Herring formula to obtain 
a completely general formula [See Eq. ( 4.38Q 1 that applies to any eigenvalue problem in the 
form of a recursion relation. This latter formula is written in terms of an action integral 
that runs between turning points, in close analogy with the continuum case. This formula 
is inconvenient for practical applications, however, and so in subsection E of Sec. |V|, we 
will transform it into another result [Eq. ( |4.39| )1 that only requires the evaluation of a small 
number of much simpler integrals. The second formula is also completely general, and has 
the advantage of making the J — > oo asymptotic structure of the splittings transparent. In 
Sec. [V|, we will apply this lattermost formula to the Fes problem, and obtain the splitting 
for all pairs of levels. We will discuss the quenching points and several other aspects of our 



results, including comparison with numerics, work by other authors |18| , ^9| , and the features 
that appear to be exact. 



II. SUMMARY OF THE DPI METHOD 

The starting point is to write Schrodinger's equation in the J z basis. Suppose is an 
eigenstate of TC with energy E. Then with J z \m) = m\m), (m\ifj) = C m , (m\TC\m) = w m , 
and (m\TC\m') = t m , m i (m ^ m'), we have 

m+2 

^ ] tm,nC- n ~\~ W m C m EC m , (2-1) 
n=m— 2 

where the prime on the sum indicates that the term n = m is to be omitted. The diagonal 
terms (w m ) arise from the Jf part of TC, the £ m , m ±i terms from the J X H X part, and the 
t m ,m±2 terms from the part. 

We can think of Eq. Q2.1J) as a tight binding model for an electron in a one-dimensional 
lattice with sites labelled by m, and slowly varying on-site energies (w m ), nearest-neighbor 
(t mj1n ±i), and next-nearest-neighbor (t m , m ±2) hopping terms. Since we can think of dynamics 



4 



in this model in terms of wavepackets, it is clear that there is a generalization of the usual 
continuum quasiclassical or phase integral method to the lattice case. This is the DPI 
method. 

Previous work with the DPI method 



11-14 has been limited to the case where the 



recursion relation has only three terms. New features arise when five or more terms are 
considered. In particular, we encounter nonclassical turning points, i.e., turning points at 
m values other than those at the limits of the classically allowed motion. It is these turning 
points that give rise to oscillatory tunnel splittings, so that this effect is absent in systems 
described by three-term recursion relations. 

The general formalism of this method [13| and the extension to five terms is discussed 
at length elsewhere ||15|| , so here we will only give a brief summary. The fundamental 
requirement for a quasiclassical approach to work is that w m and t m ^ m ± a [a = 1, 2) vary 
slowly enough with m that we can find smooth continuum approximants w(m) and t a (m), 
such that whenever m is an eigenvalue of J z , we have 



w(m) = w m , 

t a (rn^) (tm,m+a tm,m— et)/^> Ot 1, 2. 



We further demand that 



dw ( w(m 
dm 



dig _ ( t a (m,y 

J j ' dm 1 J 



(2.2) 
(2.3) 



(2.4) 



with m/J being treated as quantity of order 1. We will see that for Eq. ( [Lip , these conditions 
will hold in the semiclassical limit J ^> 1. 

Given these conditions, the basic approximation, which readers will recognize from the 
continuum case, is to write the wavefunction as a linear combination of the quasiclassical 
forms 



1 



vim) 



exp ( i / q{m')dm' 



(2.5) 



where q(m) and v(m) obey the equations 

E = w{m) + 2t\{m) cosg + 2t 2 (m) cos(2g) = H sc (q, m), (2.6) 
v(m) = dHsc/dq = —2 smq(m)(ti(m) + 4t 2 (m) cosg(m)). (2.7) 

Equations (|2.6[) and ( |2.7|) are the lattice analogs of the eikonal and transport equations. 
Equation (|2.5|) represents the first two terms in an expansion of log C m in powers of 1/ J. 

As in the continuum case, the approximate DPI wavefunction is invalid at turning points. 
These points arise whenever the velocity v (m) vanishes for given energy E, for then the 
approximation (|2.5| ) diverges. We see from Eq. (|2.7|) that v (m) can vanish either because 
sing = 0, i.e., q = or q = ir, or because q = g* = cos _1 (— ti/4t 2 ). Substituting these values 
of q in the eikonal equation, we see that a turning point is obtained whenever 



E = U (m), U w (m), or [/* 



m) 



(2.8) 



where 
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U (m) = H sc (0,m) = w(m) + 2t 1 (m) + 2t 2 (rn), (2.9) 

U v {m) = H ac (n,m) = w{m) - 2t 1 {m) + 2t 2 (m), (2.10) 

t 2 (m) 

U*{m) = H sc {q*,m) = w{m) - 2t 2 (m) - -rrV4- (2.11) 

Note that at a turning point, both m and g are determined. If we denote the values of these 
quantities generically by m c and q c , m c may be regarded as being fixed by Eq. ( |2.8| ), and q c 
by the corresponding condition q c = 0, q c = n, or q c = q*(m c ). 

To understand the nature of these turning points, let us assume that t 1 < 0, and t 2 > 0. 
[This is the case for the Hamiltonian ( |1 . 1|) . We can always arrange for t\ to be negative by 
means of the gauge transformation C m —>■ (—l) m C m . Thus there is only one other case to 
be considered, namely, t\ < 0, t 2 < 0. This is discussed in Ref. |I5| .1 It then follows that 
U n > Uq, and that 

U (m) - U*{m) = —L-faim) + U 2 {m)f > 0. (2.12) 
4i 2 (m) 

Secondly, let us think of 7i(q, m) for fixed m as an energy band curve. Then U n is always 
the upper band edge, while the lower band edge is either Uq or [/* according as whether 
—ti/4t 2 is greater than or lesser than 1. To deal with this possibility, it pays to introduce a 
dual labelling scheme for all three curves Uq, U n , and U*. We write U v (m) = U + (m), and 

= Ui(m), U*(m) = C7_(m), if q* G (0,vr), (2.13) 



Uo(m) = Ui(m), U*(m) = U-(m), if g* G (0,7 
U (m) = U-(m), U m (m) = U f (m), if g* G" (0, : 



(2.14) 



The subscripts + and — denote upper and lower band edges, while the subscripts i and / 
denote internal and forbidden respectively, since in the first case above, Uq lies inside the 
energy band, while in the second case, U* lies outside. As examples of these curves for a 
symmetric recursion relation, we show those for Fes m Fig- @- A magnified view of the lower 
left hand portion of this diagram is given in Fig. 

Turning points where E = U + , or E = [/_ when Z7_ = Uq, are analogous to those 
encountered in the continuum quasiclassical method, since the energy lies at a limit of the 
classically allowed range for the value of m in question. Points where E = U- when U- = £/* 
are physically analogous, but mathematically different since the value of q c is neither nor tt. 
Points where E = U (see the energy E\ in Fig. ^ e.g.) are novel in that the energy is inside 
the classically allowed range for m c , but the mathematical form of the connection formulas 
is identical to the case E = C/_ = Uq since q c = 0. Most interesting are the turning points 
with E = Uf (the point m = —mi in Fig. |3|, for instance), since now the energy is outside 
the allowed range for m = m c , and the value of q c is therefore necessarily complex. These 
points lie "under the barrier" and turn out to be the ones of importance for understanding 
oscillatory tunnel splittings. 

The above discussion shows that the curves Uq, U w , and U* collectively play the same 
role as the potential energy in the continuum quasiclassical method. We refer to them as 
critical curves. We have already noted that U n > Uq > U*. Let us suppose that the case of 



equality in Eq. ( 2.12 ) occurs at m — m*. Clearly t\{m*) / 4t 2 {m*) = — 1, which is precisely 
the condition found above for the lower band edge to change from q = to q = q*. Secondly, 
expanding t\ and t 2 about m*, we see that Uq and U* have a common tangent when they 
meet. 
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III. HERRING'S FORMULA FOR FIVE-TERM RECURSION RELATIONS 



The problem of computing tunnel splittings in a symmetric double well potential in 
the continuum case is greatly simplified by use of Herring's formula [ 16), PTI . An entirely 



analogous formula can be derived in the discrete case [20,14) (c) following the simplified 



treatment of Landau and Lifshitz 21 



We have already noted the importance of the critical curves. For low lying energy levels, 
in particular, the curve U- is very much like the potential energy in the continuum case, 
and it is clear that we will have an entire series of approximate energy eigenstates with 
wavefunctions localized in any one of the two wells, in the vicinity of ±m , the minima of 
U-{m). (See Fig. |.) Let C m be the nth such wavefunction localized in the right hand well, 
normalized to unit total probability, and let it satisfy Schrodinger's equation (O) with an 



energy Eq for all values of m well to the right of the left well, including in particular the 
region around m = 0. More precisely, we take C m to decay away from the right well in 
both directions. Such a function could be obtained, e.g., as the energy eigenfunction of a 
modified problem in which the on-site energy is increased by a large positive amount for all 
m < m a , where —mo <C m a <C 0, it being understood that m a is far away from all turning 
points for the energy concerned. However, this problem need not be solved explicitly, as the 
exact behavior of C m near m = —m is never needed, and therefore need not be examined 
too closely. 

Given such a function, Herring shows that the true symmetric and antisymmetric eigen- 
functions, Sj and aj, with energies E\ and E2 respectively, are given very accurately by 

f (3-D 

The product C m C_ m is exponentially small everywhere, these functions are normalized to 
unit total probaility to exponentially high accuracy. 
The Schrodinger equations obeyed by C m and 



[w m — E )C m + ^2 tm, n C n = 0, (3.2) 

n=m— 2 
m+2 

(w m -E 1 )a m + t mtn a n = 0. (3.3) 

n=m— 2 



Let us now define m r to be 1 if J is integral, and 1/2 when J is half-integral. Multiplying 
Eq. ( |3.2|) by a m , Eq. ( |3.3| ) by C m , and summing over m from m r to J, we get 

(E 1 -E ) J2 C m a m + Si - S 2 = 0, (3.4) 

m=m r 

where 

J m+2 

Si = J]] J]] dmtm,nC n , (3-5) 
m=m r n=m~2 
J m+2 

(3.6) 

m=m r n=m-2 
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To simplify Eq. 



we first note that by Eq. (|3.1| 



X! C m a m ps _ C m ^ 

m=ra r V * m=m r V * 



(3.7) 



since the product C m C_ m is everywhere exponentially small, and since is concentrated 
almost completely in the right well. Secondly, most of the terms in the sums Ex and E 2 
can be seen to be identical by shifting the summation indices in various terms suitably, and 
making use of the symmetry t m>n = t n ^ m . For example, the difference between the terms in 
Ei with n = m + 2, and those in E 2 with n = m — 2 equals 

J J J -2 

^ ] ("m^m,m+2^m+2 ^m^m,m-2"m-2) ^ ] ^m^rrt,rrt+2Cm+2 ^ ] Cm+2^ra+2,m"m 

m=m r m=m r m=m r —2 

= a m r — l^m r -l,m r +lCm r +l a m r -2^m r — 2,m r C*m r ; (3-8) 

where we have made use of the obvious facts that ij,j+ 2 and are identically zero. The 

differences between the other terms in Ex and E 2 can be similarly evaluated, and reduce to 
a small number of terms involving the product of an a with a C, which can then be written 
entirely in terms of C's using Eq. ( |3.1| ). Finally, we can see that E\ — Eq = Eq — E% = 
±A/2, and the net result is that upto an irrelevant over all sign, 



A 



to,iCo(Ci — C_x) + to,2Co(C2 — C-2) + ^-l,l(Cl 



2 t_i ifc? -C 2 +4 t_a 1 fCiCa -C_iC 

2'2 \ 9 — 9 / 2'2\2 2 2 



integer J, 
half-integer J. 



(3.9) 



Herring gives a more careful justification of his formula by employing the Temple-Kato 
error bound on energy eigenvalues |22]j23[1 . His argument can be adapted word for word to 
the present problem, and shows that the error in the splitting as calculated via Eq. ( |3.9| ) 
is exponentially smaller than the splitting itself, by a factor such as e~ cJ where c > 0. As 
J — > oo, therefore, Eq. (|3.9j) is asymptotically correct. 



We remind readers that Eq. (3.9) is not limited to the ground state splitting. 



IV. GENERAL FORMULA FOR TUNNEL SPLITTING IN TERMS OF ACTION 

INTEGRALS 

To apply Herring's formula (|3.9| ) to the Hamiltonian (|1 . 1| ) , we will use the DPI approx- 
imation for the wavefunction. Actually, we will take C m in Eq. ( |3.9| ) to be localized in the 
left well. This can only change the answer by a sign, which is not of interest to us anyway. 



A. DPI form near potential well minimum 

Let us first take up the problem of finding the DPI approximation to C m in somewhat 
general terms. Step 1 is to find C m in the classically allowed region near —mo, the minimum 
of U-(m). (See Fig. [|) We assume, as will be seen to be true for Eq. ( |1.1| ), that in this 
region U- = Uq. For energies close to U-(— m ), and m close to — m , the eikonal equation 
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can only be satisfied if q is close to zero. We can therefore expand 7i sc in powers of m + mo 
and q: 



H sc (q,m) » LL(-m ) + 



_ g 2 + _ Ma; 2 (m + mo) 2 + 



where 



M 



-[2t 1 (-m ) + 8t 2 (-m )] 1 > 0, 

-2(t! 



4U 



dm 2 



(4.1) 

(4.2) 
(4.3) 



-m 



Note that by virtue of Eq. ( |2.4| ), and its natural extension to second derivatives, u is of 
order 1/J relative to t\ and ti- 

The allowed eigenvalues and eigenfunctions can now be written down very simply by 
noting that the eikonal equation is also the Hamilton- Jacobi equation with q = dQ/dm, $ 
being the action. Thus the problem is identical to that of a harmonic oscillator. (Alterna- 
tively, we could arrive at the same result by approximating the original recurrence relation 
by a differential equation in the vicinity of —mo.) For the nth state, therefore, 



E = U-(-mo) + (n + h)u , 



and 



2 2n (n\) 2 n£ 2 



-1/4 



e- x2 ^ 2 H n (x/0, 



(4.4) 



(4.5) 



where x = m+m , H n is the nth Hermite polynomial, and £ = (Mujq)^ 1 ^ 2 . The wavefunction 
is already normalized, and the additional tails from the forbidden region only modify the 
normalization by an exponentially small amount. 

It is apparent that the expansion (|4.1| ) is invalid unless the point — m is sufficiently 
far from the edge m = — J. Since the width of the wavefunction ( |4.5| ) is y/n£, a necessary 
condition for the validity of our procedure is 



J — m > \/n£. 



(4.6) 



If this condition does not hold, then the recursion relation must be solved near the edge 
by a different method, which is tantamount to using the Holstein-Primakoff or Bogoliubov 
transformation. An example of the latter approach is given in Sec. IV of Ref. |24| . 



From the viewpoint of the DPI method, we have two turning points very close to —mo, 
one to the left, and one to the right, since the condition E = C/_(m) is then satisfied. The 
one to the left has been discussed above. Let us now consider the one to the right, and 
denote it by — m t . We have 



- m t + m 



2n + 1 



L Mu c 



1/2 



(nJ) 



1/2 



(4.7) 



The neglected terms in Eq. ( f4.1| ), on the other hand, are of relative orders g 4 , (m + m ) 3 / J 3 , 
and (m + m )q 2 / J, and thus smaller than na> for x (nJ 2 ) 1 ^ 3 . Thus, provided n <^ J, the 
solution (|4.5|) holds well past — m t , and can be matched directly onto the DPI solution under 
the barrier, without any need of connection formulas at m = —m t [25|. This argument is 



given at greater length in Sec. V of Ref. |24 . 
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B. DPI form in ordinary forbidden region 



Step 2 is to consider the DPI solution for m > —m t . Since we want this solution to decay 
as m increases, we take it as 



C„ 



B 



v(m)\ 



exp 



K(m')dm' 



(4.8) 



where n(m) = Img(m) > 0. This solution must be matched on to ( |4.5|) to determine B. We 
can continue to use the harmonic oscillator approximation (|4.1|) to 7i sc for this purpose, and 
a simple calculation [p6|, which may in fact be traced back to Furry f2"Tfl , leads to 



B 



9n 



( u g r , 



V 2tt 
n ! 



1/2 



(»+*)' 



-(n+I) 



The quantity g n is defined so that g n — > 1 as n — > oo; g = I 71 "/ 6 ) 172 
<? 2 « 1.017, .... 



(4.9) 
(4.10) 

1.075, & « 1.028, 



C. DPI form in central region 

Step 3 is to find the wavefunction in the central region near m = 0. This is already done 
if there are no turning points between —mt and m — 0. For the Hamiltonian (|1.1|) , it turns 
out that we encounter another turning point where E = Uf(m) (the only possibility) at an 
intermediate point m = —mi (see Fig. |3|). The solutions for m < —mi and m > —mi must 
therefore be related by a connection formula. To understand this turning point, we note 
that the eikonal equation fl2.6|) may be solved for cos q as 



t 2 
1 1 



, , -ti{m)±\t 2 (m) - U 2 {m)f(m]\ 1 / 2 
4t 2 (m) 

where /(m) = w(m) — 2t2(m) — E. Since cosg = —ti/4&2 at m = —mi, the discriminant in 
Eq. ( |4.11|) must vanish, and we conclude that as we cross -mi, cosg changes from real to 
complex, and q changes from imaginary to complex. [Incidentally, it may be verified that 
the condition for vanishing discriminant, i.e., 

\{m) = 4t 2 (m)(w(m) - 2t 2 (m) - E), (4.12) 

is identical to E — £7*(ra).] Since the recursion relation is real, the solution C m must also 
be real for all m. A single DPI solution for m > —mi cannot meet this demand, and so we 
must take a linear combination of two DPI forms with complex conjugate g's. If we write 
these as 

<?i,2(m) = m{m) ± y(m), (4-13) 

with k and y both real, then we must still have k > in order that C m continue decaying, 
and we may also take y > 0. Let us further write the solution (|4.8|) for m < —mi as 
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C,, 



.4 



2y/\v(m)\ 
A = 2B exp 



exp 



—mi 



—mi 
mt 



K{m')dm! 
K,(m')dm' I . 



Then, as shown in Ref. |15| , the DPI solution for m > —mi is given by 

.4 



C m — R,e 



exp (i [ qi(m')dm!\ 
Jsx{m) ^ -'~ m i ' 



where s\{m) = —iv(qi(m)) p8| . 



(4.14) 
(4.15) 



(4.16) 



D. Herring's formula with DPI approximation 



The solution ( [4.16| ) is ripe for substitution into the Herring formula ( |3.9| ). To do this, 
we first note that 



It then follows that 



cosh k cos x = — ti/4£ 2 , 
sinh/tsinx = (4t 2 / - tj) 1/2 /At 2 . 



Si = 8t2(m) sinhK(m) sinx(m) singi(m) 



(4.17) 
(4.18) 



(4.19) 



We now substitute Eqs. ( 4.16| )-( [4.19j ) into Herring's formula, Eq. ( |3.9[ ). In doing this, we 
may neglect the variation of quantities t a (m), q(m), and v(m) among the sites near the 
center of the lattice, since the number of sites involved is of order 1, and so the variation 
leads to higher order corrections in powers of 1/ J. To save writing, we denote quantities 
evaluated at m = by a bar: qi(0) = gi, k(0) = R, etc. We thus get 



e «(0+mgi) 

C m = ReA 2 , . - ■ 

V sm qi 



(4.20) 



where, 



Q — qi(m')dm', 

J —mi 

A 2 = (8t 2 sinh/€sinx) _1/2 -4. 



(4.21) 
(4.22) 



The cases of integer and half-integer J are best tackled separately. Doing the former 
first, we have 



C x - C_! = %A 2 
d + C_i = A 2 



sm qi — c.c. 



in / cos gi 
e \ . _ +c.c. 
V sm g! 



(4.23) 
(4.24) 
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C (Ci - CLi) = <^ 
C 2 - C_ 2 = 2zA 2 

Cq{C 2 — C-2, 

C 2 - C\ = iA 



e 2iQ _ e -2Imf2 / Sm Ql 

. V sin ?i 



c.c. 



zA 2 



e l!ri cos gi w sm q\ — c.c. 



e 2iQ cosg! - e - 2lmQ cosg^ 



/ sm gj 
sin qi 



cos ft e^cosgx-e- 21 ^,/^^ 
V V sm 9i 



c.c. 



c.c. 



Substituting these and the formula t\ = — 4t 2 cosh k cos \ m to Eq. (|3.9|), we get 



-2Imfi sm Ql 



sm qi 



Re9 



A = -8A 2 t 2 Im \^e 2lU 6 - < 

= cos ft — cosh R cos 

But, it follows from Eq. fl£T5|) that 

cos ft = cosh k cos x — i sinh R sin x, 

so the second term in Eq. (|4.29|) vanishes altogether, and 

A = AA 2 2 t 2 sinh R sin x(e 2in + e~ 2in * 
= ±A 2 {e 2in + e~ 2in *), 

where we have used Eq. ( |4.22|) in the last step. 
For half-integer J, we get 



A 2 

n 2 A 2 

A 2 

C±l/ 2 C±3/ 2 = — 



_2ifi „-2Imn 



+ c.c. 



^ sm qi 

^2in 



sm q 1 1 ) 

-2Imn 



f " e ±2l?1 + e^" 1 I + c.c. 



sm qi 



sm g! | 



(4.25) 
(4.26) 
(4.27) 

(4.28) 



(4.29) 
(4.30) 

(4.31) 



(4.32) 



(4.33) 
(4.34) 



Thus, 



and 



^2 ^2 1 a2iJ1KI 

u l/2 ~~ °-l/2 — 2 2 



C 1/2 C 3/2 - C_ 1/2 C_ 3/2 = L4 2 (cosgie 2iQ - c.c. 



A = iA\ 



(fx + 4t 2 )e 2in - c.c. 



AA\t 2 sinh « sin X {e 2in + e~ 2iQ * 



(4.35) 
(4.36) 



(4.37) 



which leads, once again, to Eq. ( |4.32j ). 

Collecting Eqs. ( |4.9| ), ( |4.15| ), ( |4.21|) and ( |4.32|) , and making use of the symmetry of the 
problem, we may write the tunnel splitting for both integer and half-integer J as 
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A 



2tt 



exp 



q(m')dm! ) + c.c 



■mt 



(4.38) 



Here q(m') is chosen to have a positive real part x i n the first term. We note once again 
that this result applies to higher pairs of excited states, and not just the ground pair. The 
essential dependence on n, the excitation number, enters through the n dependence of m t , 
the turning point. 

The similarity of Eq. ([4.38 ) to the final result in Ref. is striking ||29|| , and one can 
ask whether one should not have anticipated it right away. For the ground state pair, the 
instanton approach [0(a) makes it very easy to understand the presence of two complex 
conjugate tunneling actions, and the fact that they should be superposed, but does not 
give the prefactor. The action integrals in the instanton approach, however, run not from 
turning point to turning point but from one minimum of the energy to the other. Further, 
properly justifiying the prefactor using instantons has proven very difficult 1 30 ,|31| . Purely 
as a recipe for calculations, however, a hybrid approach, in which one adds the tunneling 
actions from all equivalent instantons, and uses the DPI approach to determine the form 
of the prefactor, would appear to be valid for all problems. Thus, we strongly suspect that 
Eq. ( f4.38p is correct even when the recursion relation has seven or more terms. 



E. Extraction of singular parts of action integrals 

While the formula fl4.38|) is very general, it has the disadvantage that the action integral 
runs between turning points. The integrand is therefore close to a singularity, and for low 
lying states, this gives rise to terms in the action that depend on In J. Hence the formula 
does not reveal the asymptotic behavior as a function of J in a transparent way. 

In this subsection, we will show that we can write the splitting very simply in a way that 
does not suffer from the above drawback. The final result is 

A n = lW-^ F" + 5 e - r °cosA n , (4.39) 
nl V 7T 

where, 

r° 

T = 2 K (m)dm, (4.40) 

J -mo 

A n = 2 f ( X o + (n + |VoXoV m ' ( 441 ) 

J— mi \ J 

F = 2Mto (m - mi) 2 exp ^-2 (Qi + co J° K dm^ , (4.42) 
(h=\ \^£L + ^—)dm. (4.43) 

In Eqs. ( |4.40H4~43| ) , the irregular turning points imi may be evaluated by setting E = 
U-(±m ), and it should be recalled that ±m are the minima of U-{m). Further, 
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k, = K(m, e = 0); k' c 
Xo = x(m,e = 0); x'o 



dK(m, e) 



de 
dx(m, e] 



e=0 



e=0 



Bq = cos q(m, e = 0); B' 



de 

d cos g(m, e) 



5e 



e=0 



with 



e = E — UJ-mo). 



(AAA) 
(4.45) 
(4.46) 

(4.47) 



The problem of finding the low level splittings is thus reduced to the evaluation of 
a handful of integrals. The proliferation of notation masks the actual simplicty of these 
formulas. 

To derive these results, we follow a procedure similar to that used for the continuum 
case in Ref. ||26|| . We begin by defining 



$(e) 



■mt(e) 



q(m, e)dm, 



(4.48) 



where the energy dependence is made explicit. The splitting for the nth pair of states is 
then given by 



^( e -»fcO +c .c.), 

Z7T 



(4.49) 



with e n = (n + \)uJq. Writing x = m + m as in Eq. (|4.5|), the integrand in $ behaves as 
(x 2 — a^) 1//2 near the lower limit, with x t = — m t +m ~ e 1 / 2 . Thus there is a singular part in 
$ of the form elne, which it is our goal to extract. To this end, we differentiate Eq. ([4.48 ) 
to get 



3 <9g 

—am. 

m t (e) de 



(4.50) 



Note that the term arising from differentiating the lower limit vanishes, nor is there any 
explicit contribution from the singular behavior q ~ (m + m^) x l 2 for m near — m c . 

Next, let us divide $'(e) into two integrals, in which the limits of integration are 
— m t and —mi, and $' 2 , which runs from —mi to 0. Defining 



B e (m) = cos(g(m, e)), 



we have 



=dm, 



-"h(«) ^B 2 e {m) - 1 
where = dB e /de. It follows from Eq. ( |4.1| ) that near m 

5 e w 1 + (|Mcj 2 x 2 - e) M + • • • . 
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(4.51) 



(4.52) 



(4.53) 



so the integrand in Eq. ( |4.52|) behaves as — \/ojq{x 2 — x 2 ) 1 ^ 2 . If we add and subtract the 
integral of this expression, we obtain 



j-i 



dx 



+ 



x+ 



Xi 



B> 



B 2 (m) - 1 



+ 



U}Q\ x 



.1 4 



(4.54) 



where x\ = m — rri\. The first integral can be evaluated exactly. In the second integral we 
can put e = both in the limits and in the integrand, since we are not interested in terms 
of 0(e). Ignoring terms of this order throughout, and making use of the relation 



.1 I. 



2e/Mcu t 



we obtain 



2lJn 



hi 



2Mu (mo — mi) 2 



+ 2Q, 



(4.55) 



(4.56) 



where Qi is given by Eq. (|4.43 ). Also, we can evaluate mi at e = 0. 

The remaining contribution to &(e), ^(e), can be evaluated simply by putting e = 0, 
since the neglected part is 0(e). Recalling the definitions ( |4.44| ) and ( [4.45[ ), we have 



(4.57) 



Oil 



We now integrate the expression for $'( e ) an d obtain $. It is useful to separate the real 
and imaginary parts of the answer at this stage. For the real part, we get 



2 Re$ = T + — 

U 



2Qi - 1 + In 



2MLU (m — mi) 



+ 2Ur 



k'q dm 



(4.58) 



with To given by Eq. (|4.40|) , while for the imaginary part, A n = —2 Im$, we get Eq. ( [4.41 ). 

Substituting Eqs. ( |4.58| )-( f4.41| ), and the definition ( |4.10| ) of g n in the formula ( [4.49| ) for 
A n , and recalling that e n = (n + |)co>o> we finally obtain the answer quoted at the start, 

Eq. flogp . 



V. APPLICATION TO THE Fe 8 PROBLEM 

We now apply our general formulas to the specific problem of Fes, as described by the 
Hamiltonian ( |1 . 1|) . The various matrix elements of this Hamiltonian are given by 

w m = ^(h + k 2 )[J{J + 1) - m 2 ], (5.1) 

t m , m +i = -^gfi B H x [J{J + 1) - m(m + 1)] 1/2 , (5.2) 

t m ,m+2 = \{h - k 2 ) [[J(J + 1) - m(m + 1)][J(J + 1) - (m + l)(m + 2)]] 1/2 . (5.3) 

We must now replace these by continuous functions w(m), ti(m), and t2(m). Since our 
formalism requires knowing the first two terms in the action in an expansion in powers of 
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1/J, it follows that we need only determine the functions w(m) etc. to the same order. 
Furthermore, this determination need not be made in the form of a power series, and any 
functional representation that gives the first two terms correctly will be adequate. The most 
convenient way to do this is to replace the combination J( J + 1) in the above expressions 
by J 2 , where 



J = J 



1 

2' 



(5.4) 



The evaluation of the integrals fl4.50|) - (|4.43 ) is then lengthy, but straightforward. We will 
present and discuss the final results first, and give the details of the analysis later. 



A. Tunnel splittings for Fe§ 

The final result for the splitting of the nth pair of levels is 



1 



A n = -W-a; F n +5 e -*o cosAn) 
n! V 7r 



where, 



u 
F 



2J[hk 2 (l - /4)] 1 / 2 , 

A^fl- hi) 3 / 2 



8J- 



1 - A - hi 



J 



In 



'1 - hi + vT 



'l-hl-^/X, 



dn 



[1 - hl)(l - X) + h x V\" 



\J(l-hl)(l-\)-h x V\, 



A, 

In Eqs. 



max < 0, 7T J 1 



VT^XHJ 

X = k 2 /ki, and 

J H T 

h v 



mx 



h x o 



H x 



(5.5) 

(5.6) 
(5.7) 

(5.8) 

(5.9) 



(5.10) 



Recall that H c = 2k\J / ' qhb- 

Let us now turn to the discussion of these results. The first point concerns the fields 
where the nth tunnel splitting vanishes. Taking account of the fact that A n is necessarily 
positive as indicated by Eq. (|5.9| ), we see that this happens whenever ||[!|] (d,e) 



Hx 



VT^X 
J 



J 



(5.11) 



with I = n, n + 1, 2 J — n — 1, yielding 2(J — n) quenching points in all for A n . When 
n = 0, these are the results quoted in Sec. |. 

In Fig. |4] we compare Eq. ( |5.5| ) with the numerically evaluated splittings for the first 
three pairs of levels. Within our numerical precision, we always find the zeros of A n to agree 
with Eq. (|5.11|) . Note, however, that for other values of H c , the discrepancy between the 
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numerics and Eq. ( |5.5|) is well outside our numerical error, so that Eq. (|5.5|) is not exact, 
even though as an asymtotic estimate of the splitting it is rather good. This means that the 
leading semiclassical approximation does not give the eigenvalues themselves exactly, and 
only the quenching points appear to be so reproduced. The second point to note is that for 
n = 1 (the pair of first excited states in each well), the highest field quenching point is lost, 
for n = 2, the highest two points are lost, and so on, exactly as indicated by Eq. ( |5.11| ). 

Next, let us compare our answers with previous work. Let us consider the Gamow factor 
r first. Except for the replacement of J by J and H x q by h x , this is precisely the action in 
Eq. (3.10) of Ref. [|T|(c). This agreement is unsurprising, because if we write A n in the form 
of a prefactor c\ times a Gamow factor exp(— Jc 2 ) where c 2 = 0(1), then the J —>■ J, and 
h x o h x corrections in Eq. (|5.8|) represent terms that should be included in the prefactor 
Ci, which we did not seek to find in Ref. 0](c). The detailed form of the prefactor is perhaps 
more interesting. Up to multiplicative terms of order J°, our answer for A n agrees precisely 
with that in Ref. |TB],[]1| . We do not understand, however, how these papers have succeeded 
in sidestepping the difficulties in the path integral treatment that were noted by Enz and 
Schilling [ftC|| , and by Belinicher, Providencia, and da Providencia J3T]. In Ref. [[I1J, for 



instance, the problem is treated by writing the spin coherent state expectation value of the 
Hamiltonian ( |1.1D in spherical polar coordinates, and integrating out cos 9 (the J z projection) 
exactly, and then addressing the resulting effective Lagrangian for the <fi coordinate exactly 
as for a massive particle in one dimension. In performing the integration over 9, however, it 
is not clear to us why S 2 is replaced by S(S + 1) in the scalar potential V(<p) [see Eq. (12) 
there], but not in the vector potential 0(0). 

A related point, which is relatively minor, but has scope for creating confusion, is that 
if the Gamow factor is written as exp(— Jc 2 ) with c 2 = 0(1), then it is safest to write the 
J dependence of the prefactor as cu J n+1 ^ 2 , since u depends on parameters such as k± and 
fc 2 , whose scaling with J is a matter of choice, at least as far as model Hamiltonians are 
concerned. 

One further check is obtained by considering the limiting case H x = 0, answers for which 
are known. [See, e.g., Eq. (16) of p0| , Eq. (48) of or [^4j].] Transcribing Eqs. (4.30) and 



(4.31) from MM in terms of the present parameters, we get 



A n = -^F "A ; (5.12) 
F Q = SJj^v ( 5 - 13 ) 



jxi/2 A i/4 /i_VA 



7T 



Ao = 8w M - 7= -7=) , (5.14) 



with ujq = 2J(/ci£; 2 ) 1 / 2 . It follows from Eqs. (|5.6|- |5.9|) that as H x — > 0, cu — >• cu 00 , F — > F 
[see Eq. ( |5.7|) 1, cosA n — > ±1, and 

J + !)l„i±^. (5.15) 



It is then easy to see that our present answers for A n go over precisely into Eqs. ( |5.12| )-( |5.14| ). 
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B. Evaluation of Action Integrals 



To carry out the evaluation of Eqs. ( |4.5U )-( 4.43[ ), it is convenient to measure energies 
(including u ) in units of fciJ 2 , and introduce the scaled variable /x = raj J. In terms of 
these variables, we have 



M m ) = (l + A)(l-/i 2 )/2, 
t l {m) = -h x {l-^) l l\ 

t 2 (m) = (l-A)(l-/i 2 )/4. 



(5.16) 
(5.17) 
(5.18) 



The turning points /Jq = mo/ J, and \i\ = mi/ J (for e = 0) are given by 

0o = (1 - hlf'\ 

*ii = [(1- A- /£)/(! -A)] 1 / 2 . 



(5.19) 
(5.20) 



It is most convenient to express everything in terms of //q and fi±, so we give inverse formulas 
as well: 



A = (0o-0i)/(l-0?)- 



(5.21) 
(5.22) 



The mass and the small oscillation frequency are given by 



M 



2\hl 



u = 2[\(l-hl)) 1 / 2 /J 



2(l-/i o 2 )(^-0? 
200 / 0o -0i * 



2\ ' 



(5.23) 
(5.24) 



To evaluate the integrals, we need expressions for k , Xo? ^o> e ^ c - m the ranges 0i < 
/i < fj,Q, and < /x < 0i- The requisite calculations are straightforward so we give the main 
results only. First, in the range H\ < n < /x , we get 



Bq = cosh Kq 



2 -[(^-0p(0 2 -0D] 1/2 
[(i-^)(i-0 2 )] 1/2 



5 2 -i 



^0 

5o 2 -l 



sinh k 



1 



00 — 01 



2 -0i 



0i 



[(l-0g)(l-0 2 )] 1/2 



1-0? 



2[(l- A , 2 )(l- / i2)( /i 2_ /i 2 )(;u2 _ /i 2 )] l /2 : 
00 1 



J 



2 - 01 



0o-0i 



2 -01 



(5.25) 

(5.26) 
(5.27) 
(5.28) 



Next, in the range, < \x < 0i, we first put e = in Eqs. ( [4.1 7[ ) and (|4.18[ ), and solve, to 
obtain 
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cosh« = [(l-/iMl-^o)] 1/2 , 

cos X o=[(l-^)/(l-^ 2 )] 1/2 - 
sin Xo = [(^-/i 2 )/(l- / i 2 )] 1 / 2 

We then differentiate Eqs. (|4.17|) and (f4.18|) , and set e = to obtain the equations 

sinh k cos Xo K o ~ cosh k sin xo Xo = 0> 

cosh n sin Xo + smn K o cos Xo Xo = ~ V 8*2 sinh k sin xo- 

Solving these, we obtain 



Xo 



2(^-/x 2 ) v> 2 ~/i 2 )- 1/2 



(5.29) 
(5.30) 
(5.31) 



(5.32) 
(5.33) 



(5.34) 



The first integral that we wish to evaluate is To, which will give us the dominant WKB 
or Gamow factor in the tunnel splitting. We divide the integral into two parts by breaking 
the integration range at mi. From the right-hand part, an integration by parts gives 



r i = 2j 

= 2 J 



/'() 



Kodjj, 



K o(/^C - 



w fi dB (fi) 
4 sinh k dfi 



dfi 



(5.35) 



while from the left-hand part we get 

r 02 = 2J K dfi = 2Jk (/j, 1 )hi, (5.36) 
Jo 

as k,q is a constant in this range. Since ^0(^0) = 0, To2 cancels the first term in Eq. ( p.35| ), 
leaving us only with the second for Tq. Using Eqs. ( |5.25|) and (|5.26| ), we find 



T Q = 2J{l-fi{ 



2x1/2 



dfi 



mi 



/i 2 )(^ 2 -/ii) 1/2 ' 



(5.37) 



The integration is now elementary, and the result, expressed back in terms of A and h x is 
Eq. Q. 

The second integral to be evaluated is A n . For the first term in Eq. ( |4.41| ), we integrate 
by parts, and use Eqs. fl5.30p and ( |5.31| ): 

/ Xo{^)dfi = 2 J I — — cosxo dfi 

Jo Jo 



sin xo dfi 



2J 



-dfi 



l0 (l-/i 2 )(/if-/X 2 )V2 

= 7rJ[l-(l- A i 2 ) 1 / 2 ]. (5.38) 
For the second term in Eq. ( f4.41p , we have, with e = {n + |)co> , and Eqs. ( |5.24| ) and ( p.34[ ), 



2eJ 



Ml 



n' dfi = -{2n + l)fi Q (fil - fi\) l/2 
= - (n + |) 7T. 



Mi 



dfi 



o {fil-fi>){nl-f£)V* 



(5.39) 
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Adding together the parts, and rewriting the result in terms of H x and A, we get Eq. Q5.9|) . 
The restriction that A n be positive follows from the fact that we chose q(m) to have a 
positive real part in Eq. ( |4.38| ) . Thus A is necessarily positive as defined in Eq. ( [4.41| ). If 
H x is so large as to yield a negative value for the function of H x that results after doing the 
integral, that means that in fact there are no irregular turning points in the problem. Both 
terms in Eq. (|4.38| ) are then equal, and the formula reduces to the expected one when there 
are only regular turning points. 

Note that unlike Eq. ( |5.8| ), what appears in Eq. ( |5l)|) is the ratio H x /H c , i.e., h x0 , not 
h x . This fact is important for the location of the diabolical points. 

The third integral we need is that of n' from to H\. Using Eqs. ( |5.24| ) and ( |5.34| ), we 

get 



- /-mi /-mi dfi 

2u J / K dfi = -2/i / -5 5- 

Jo Jo fiQ — m 



In 



Mo ~ Mi 
/i + Mi 



(5.40) 



The fourth and last integral needed is Q\. Substituting Eq. ( |5.28| ) in Eq. ( f4.43| ), we 
obtain 



Qi 



mi 

Mo 



Mo 



M 



Mo — Ml 



M 



Mi 



Mo - M 



dfi. 



(5.41) 



The integrand is now nonsingular at fi = fi . We can make this manifest by rationalizing 
the difference of square roots in the first term. Some simple algebra yields 



Qi = - 



Ml 



Mo 



Mo + M 2 - M? 



/i 2 - [i\ no 



Mo 



Mi + MvM 



Mi 



We now make the substitution \i = \i\ cosh z, and define 

cosh^o = Mo/Mi- 



(5.42) 



(5.43) 



This yields 



2 o cosh 2zq + cosh 2z 



sinh 2zo + sinh 2z 
z o cosh(z + Zq 
sinh(z + 2 



= — ln(2 cosh z ) 



-In 







2/i, 
Mi 



(5.44) 



We now have all the ingredients needed to calculate the quantity F. Substituting 
Eqs. (g]23|), (|5T2lD , QCTD , and ( gg§ in Eq. (^742|) , and writing the result in terms of A 
and /ia,, we obtain Eq. (|5.7|) . Note that in writing down the final answer, we have replaced J 
by J and h x by h X Q in this formula. This is because F is part of the pre-exponential factor 
in A n , which is determined only to leading order in 1/J. Keeping higher order corrections 
by distinguishing between J and J or h x and h x0 is not justified. 

The final answer ( |5.5| ) for A n is obtained by substituting Eqs. (|5.8|), (|5.9|), and ( |5.7| ) in 
Eq. ( |Q9D . 
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FIGURES 



FIG. 1. Spectrum of the Hamiltonian (1.1) for J = 3, as a function of H x /H c . H z /H c = 0, 
0.07454, and 0.1491 in (a), (b), and (c) respectively. The small ovals indicate points that are 
narrowly avoided anticrossings, but appear to be crossings on low resolution. 



FIG. 2. Critical energy curves for the Hamiltonian (1.1), showing the dual labelling scheme. 



FIG. 3. Magnified view of the lower left hand region of Fig. Q showing the point of tangency 
m* between Uq and [/*, and turning points at m = —m% and —rrt\ for an energy E. 



FIG. 4. Comparison between numerical (solid lines) and analytic [Eqs. ( |5.5| -0.9|), dashed lines] 
results for the splitting between the first three pairs of levels for H z = 0. The parameters are 
fci = 0.321 K, k 2 = 0.229 K, close to those for Fe 8 . 
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